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The strong-coupling regime of Kardar-Parisi-Zhang sur- 
face growth driven by short-ranged noise has an upper crit- 
ical dimension d> less or equal to four (where the dynamic 
exponent z takes the value z{dy) ~ 2). To derive this, we 
use the mapping onto directed polymers with quenched dis- 
order. Two such polymers coupled by a small contact attrac- 
tion of strength u are shown to form a bound state at all 
temperatures /3~^ < /3^^ , the roughening temperature of a 
single polymer. Comparing the singularities of the localiza- 
tion length ~ M^'''^"''' at /J"^ and ^± ~ u^/^'^"^) below 
/3~^ then yields d> < 4. 

PACS numbers; 5.70.Ln, 64.60.Ht, 68. 35. Fx 



The Kardar-Parisi-Zhang equation has been intro- 
duced as the simplest nonlinear evolution equation 

dth{r, t) = v V'h{r, t) + ^ (V/i(r, t))^ -f r;(r, t) (1) 

for a continuous "height field" h{r, t) driven by Gaus- 
sian noise ri{r, t) with {ri{r, t)) = and {ri{r, t)ri{r' , t')) ~ 
a'^6'^{r — r')S{t — t'). It appears ubiquitously in a num- 
ber of noncquilibrium statistical problems such as fluid 
dynamics or dissipativc transport, as well as in systems 
with quenched disorder, for example flux lines in a dirty 
type-II superconductor (see Q for reviews). 

The morphology of a growing surface governed by (|l|) 
is well known in low spatial dimensions. The asymptotic 
scaling of the spatio-temporal height correlations. 



{{h{v,M)-h{v2M)f)^r^^C{tlr^) 
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(with r = |ri — and t = \ti ~ t2\) deflnes the rough- 
ness exponent x and the dynamic exponent z. In the 
Gaussian dynamics (A = 0), one has x— (2 — (i)/2 and 
z = 2. For A ^ 0, Galilei invariance imposes the re- 
lation X + z = 2 In dimensions d < 2, any small 
nonlincarity {\/2){'\Ih)'^ is a relevant perturbation of the 
Gaussian theory and leads to new values of the exponents 
(X = 1/2, z 3/2 in d = 1 g and X « 0.386, z k. T612 
in d = 2 1^). In the renormalization group, there is 
a crossover between the Gaussian fixed point, which is 
(infrared-)unstable, and the strong- coupling fixed point, 
which is stable. For d > 2, a Gaussian surface is 
smooth. A small nonlinearity (A/2)(V/i)^ does not al- 
ter this asymptotic scaling. There is now a roughening 
transition to the strong coupling regime at finite values 



±Ac iQ-j^]- In the renormalization group, the transition 
is represented by a third fixed point. This critical fixed 
point is unstable and appears between the Gaussian fixed 
point and the strong-coupling fixed point which are now 
both stable HJl^. The critical point is characterized by 
the dimension-independent critical exponents x* — and 
z* = 2|ll|,@. 

It is a notorious difficulty of the strong-coupling regime 
that above dimension one its properties are inaccessible 
to any known systematic approximation, let alone to an 
exact solution. In particular, renormalized perturbation 
theory with an e-expansion about d = 2 fails to produce a 
strong-coupling fixed point jlj]. This feature, which the 
Kardar-Parisi-Zhang equation shares with other driven 
systems such as fully developed turbulence, is a ma- 
jor open problem in nonequilibrium statistical physics. 
In this letter, we discuss a (modest) step towards its 
eventual solution. We find that the Kardar-Parisi-Zhang 
equation has an upper critical dimension d^ < 4, where 
the exponents in the strong-coupling regime take the val- 
ues z = 2 and x = 0, equal to those governing the limit 
d ^ oo |jl^,|l^ (and thus presumably the whole interval 
d> < d < oo). As d> is approached from below, the 
exponents tend to these values continuously. Hence, the 
upper critical dimension could serve as the starting point 
for a controlled expansion. However, the reader should 
be cautioned that the name "upper critical dimension" 
may be misleading since d> does not mark the borderline 
to simple mean-field behavior as in the standard theory 
of critical phenomena, but instead to an even more com- 
plicated state in high dimensions with presumed glassy 
characteristics |p^ , p^ . 

Even the existence of a finite upper critical dimension 
has been very controversial. Numerical work seems to in- 
dicate that a strong coupling phase with nontrivial expo- 
nents z < 2, X > persists in dimensions d = 3, 4, . . . , 7 
Various theoretical arguments, on the other hand, 
favor the existence of a finite upper critical dimension 
d> . This is the case for functional renormalization group 
calculatio ns pq ], which are supported by a Flory-type 
argument ]17|, and for a 1/ d-expansion based on a solu- 
tion on the Cayley-tree (criticized, however, in [^). 
Recent work treating the Kardar-Parisi-Zhang dynamics 
in mode-coupling approximation gives contradictory re- 
sults, predicting values of d> between 3 and 4 |ljjl^ or 
d> = oo 

In contrast to all these approaches, the arguments of 
this Letter are not based on any uncontrolled approxi- 
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mation scheme. We use the known equivalence of the 
Kardar-Parisi-Zhang equation to a system of directed 
polymers and compare exact (and individually testable) 
physical properties of this system in the strong-coupling 
phase and at the roughening transition. The importance 
of d = 4 for the roughening critical point has already 
been stressed recently | p^ , pO[ . For 2 < d < 4, this uni- 
versality class is accessible by renormalized perturbation 
theory, which produces the exponents x* = and z* = 2 
exactly to all orders. At d = 4, however, there are singu- 
larities in some observables (for example, in the exponent 
y* defined below). 

Via the well-known Hopf-Cole transformation. 



exp 



2^ 



^(r/,t/) 



with the Hamiltonian 



n = 









f ' dt 


\ 




'o 




\At) 



VrS{ritf)-rf)exp[-pn] (3) 
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the Kardar-Parisi-Zhang equation can be mapped onto 
the equilibrium problem of a directed polymer r{t) living 
in the quenched random potential Ary(r, t) at temperature 
f3~^ = 2v. The polymer is characterized by its transver- 
sal displacement 
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Mh)-vit,W)^\h-h\^'^ 

and by free energy quantities like the "Casimir" term 
7^(i?) = lim dT(F{T, R) ~ F{T, oo)) - r('^-'^)/( (6) 

T — >oo 

in a system of longitudinal size T and transversal size 
R. (Averages over the disorder are denoted by overbars, 
thermal averages by brackets (...).) The asymptotic scal- 
ing in (||) and (|^) is related to the growth exponents by 
^ = 1/z and uj = x/z; the scaling relation due to Galilei 
invariance now reads a; = 2^ — 1. In the strong-coupling 
regime for d < ci>, the polymer becomes superdiffusive 
(C > 1/2), and its free energy acquires an anomalous 
dimension lu > 0. These disorder-induced fluctuations 
persist in the zero-temperature limit, that is, in the en- 
semble of minimum energy paths ro(t). In the weak- 
coupling (high- temperature) regime for d > 2, thermal 
fluctuations dominate (C = 1/2) and hyperscaling is pre- 
served (w = 0). The roughening transition between these 
two phases takes place at a finite temperature /?~^. At 
d = d> (and probably for d > d>), the exponents C = 1/2 
and a; = govern the low-temperature phase as well. 

In its mapping to directed polymers, it is a natural 
extension of the model to study an ensemble of lines (t) 
that live in the same random potential and are coupled 
by direct mutual forces p^-|2^ . The simplest case is that 
of just two lines with the interaction 



Hi = -uJ dt^{t) 



in terms of the local pair field ^{t) = S{ri{t) — r2(t)). 
For d > 2, this interaction is irrelevant at the Gaussian 
fixed point governing the weak-coupling regime. How- 
ever, for d < d>, it turns out to be relevant both at 
the strong-coupling fixed point ||2^ and at the transition 
fixed point. Thus for all temperatures /3^^ < /3~^, an 
attractive contact interaction (u > 0) leads to a bound 
state, suppressing the thermally activated relative fluctu- 
ations of the lines on scales |ri — r2 1 ^ C-L- The transversal 
localization length has the singularities 



1 



- UJ 
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S,± ^ u with y 

in the low-temperature phase (/3~^ < f3~^), and 

u~^/^y' with y* ^{d- 2)/2 . (9) 

at the critical temperature ■ Of course, for any u > 
kept fixed, S^±_ should increase with temperature, since 
the amount of entropic fluctuations increases. This im- 
plies y/C > 2y* , that is, an upper bound on the free 
energy exponent, w < (4 — d)/d. Our main result d> < 4 
then follows immediately. We now discuss the arguments 
leading to Eqns. (||) and (||) above. 

(a) Strong coupling phase: At low temperatures and 
for a; > 0, two noninteracting lines {u = 0) in the 
same random potential have a stationary pair distribu- 
tion P{r) = Yant.T-t^oa {5{Yi{t) - V2{t) - r)) given by 



P{t) 



with 9 



(10) 
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in a suitable normalization |gl|,Q. The nonintegrable 
short-distance singularity is cut off on scales r ^ , 
where fx ~ Pr'^^ is the crossover length to the asymp- 
totic strong-coupling behavior [ p5[ . Thus with finite 
probability, the two lines share a common "tube" of 
width along the minimum energy path ro(t) for a 
single line; this probability approaches one in the limit 
[3^^ — > 0. In that sense, the ground state r'o(i) is 
unique At finite temperatures, however, the lines 

make thermally activated individual excursions on all 
length scales from the tube tq (t) . These excursions gen- 
erate the power law distribution ( |l0|) , all positive integer 
moments of which diverge. 

Additional interactions between the lines probe the 
uniqueness of the ground state as well as the statistics 
of the low energy excursions pst . For example, a mu- 
tual repulsion (^ with u < forces one of the lines onto 
a distant excited path. (In superconductors, this is im- 
portant to stabilize a dilute ensemble of flux lines in the 
strong-coupling regime against collapse.) For the pur- 
pose of this Letter, it is more useful to study a weak 
attractive interaction (u > 0), which localizes the two 
lines to each other. The normalized bound state distri- 
bution has the form P(r,^_L) = N-'^{S,i_)P{r)T{r 
on scales r ^ , it falls off exponentially as given by the 
scaling function T, but in the scaling regime r <C its 
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dependence on ^± originates only from the overall nor- 
malization N{^±) = J d'^rP{r)T{rl£_i_). 

Pair interactions have been treated in ref. by 
renormalized perturbation theory for the Hamiltonian 
based on the short-distance expansion '^{t)'^(t') ^ 
/3^^|i — t'\^'^^{t) at the strong-coupling fixed point of 
noninteracting lines {u — 0). For u > 0, one finds a 
bound state with the localization length singularity (||). 
This is in agreement with results in = 1 from numerical 
work 



12 23 1 and from the dynamic renormalization of an 
extended Kardar-Parisi-Zhang equation p9{ . 

It is instructive to cast the field-theoretic derivation of 
ref. |2^ ] into the form of a variational scaling argument: 
a weak bound state has a localization length that min- 
imizes its free energy per unit of t, balancing the free 
energy gain 5f^ from the overlap of the two lines with 
the loss 5/^due to the confinement of their relative fluc- 
tuations [gO|. The overlap free energy is proportional to 
the stationary expectation value of the pair field, 



5J, = ^u{^) ^-uN-\U) 



(11) 



The confinement energy has to vanish in the limit 
0, where all relative fluctuations of the lines are sup- 
pressed even without attractive forces. Therefore, in- 
stead of the temperature-independent Casimir energy 
^("-i)/C analogous to (||), one expects 



6h 



(12) 



a term that is analytic in the scaling variable (3^^ and 
respects hyperscaling |32|. The variation of + 5f ^ 
with respect to then leads to (^). 

(b) Roughening transition: In contrast to the strong- 
coupling fixed point, the renormalization group for the 
critical fixed point in 2 < d < 4 is well understood. 
The polymer partition function (H) has been shown to 
be one-loop renormalizahle |l^j20|] . It is this property 
that produces the exact dimension-independent expo- 
nents = 2 and z* = 2; these values agree with previ- 
ous results from numerical work in d — 3 p3| , from dy- 
namic renormalization group calculations to one-loo p P ] 
and two-loop [p^ order, and from scaling arguments [[ll[ . 
At the critical point Ac, small variations of A are a rel- 
evant perturbation of dimension y* — [d — 2)/2 (with 
t as the basic scale) The local field conjugate to 

A, $,,(i) EE / d'^r' ri{Y\t)5'^{Y{t) - r'), encodes the ran- 
dom potential evaluated along the polymer paths. This 
field generates the crossover from the critical fixed point 
to the strong-coupling fixed point (A ^ oo) and to the 
Gaussian fixed point (A = 0). It is now a purely formal 
matter to relate its disorder-averaged correlation func- 
tions to those of the pair field for u = p^ , p6[ . 
The simplest example is the stationary one-point func- 
tion (<i>,,)(i?) obtained from the free energy per unit of t, 



f[R) = -/?-ilimT^ooi9TlogTrexp(-/3H)(T,i?). With 
the disorder correlation and the pair interaction regu- 
larized on the microscopic scale a, rj{v' ,t')ri{v" ,t") — 
aHa{v' - v")5{t' - t") and^(t)^(ri(i)_-r2(i)), the 
universal parts [Q of ($,,(f)) = ($^) and are seen 
to be proportional: 



-A-a^|dVdV'5„(r'-r")^ 



t)Srj{r",t) 



(3a^ / d''r'd''r",5a(r' - r"){5{r{t) - r')S{r{t) - r"))' 



= -l3<^^J d^r'dV,5a(r' - r"){S{r{t) - r')){S{r{t) - r")) 

= -l3a^l^. (13) 

At A = Ac, the correlation functions of <I>^ are scale- 
invariant with the exponent x* = 1 — y* — {A — d)/2. By 
(p^), the same holds for the correlation functions of 
for example. 



(*,)(i?,Ac) ^ {%){R, Ac) 



(14) 



Hence, the local pair interaction ^a{t) is like ^nit) a rel- 
evant scaling field of dimension x* at the critical fixed 
point. This result is consistent with the one- loop dy- 
namic renormalization group discussed in ref. We 
conclude that the bound state generated by this interac- 
tion has the localization length singularity (||) at 
Just below that temperature, there is a crossover from (S) 
to (^) at the scale w ~ - /3^^ 

The singularity (|^) can again be obtained from a vari- 
ational scaling argument. At /3^^, the competing free 
energy contributions analogous to (|ll|) and ( p^ ) are the 
overlap term 



-x*/C 



(15) 



following from ( p_4D and a Casimir term of the form 

- d""'^^'* = ei^ 



(16) 



It is tempting to speculate about the nature of the 
strong-coupling regime in high dimensions. Below d>, 
the uniqueness of the ground state implies a collapse of 
the pair distibution function P(r) in the limit Pn — > 0, 
see Eq. (p^. This property may be lost in a glassy 
phase. If many states remain accessible to the system, 
long-ranged pair correlations P{r) ~ may persist even 
as f3R — > 0. Thus certain characteristics of the strong- 
coupling regime may qualitatively resemble those at the 
roughening transition below d> . 
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